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ABSTRACT 

CIZA J2242.8+5301 (z = 0.188, nicknamed ‘Sausage’) is an extremely massive (M 200 ^ 
2.0 X IO^^Mq), merging cluster with shock waves towards its outskirts, which was found 
to host numerous emission-line galaxies. We performed extremely deep Westerbork Synthe¬ 
sis Radio Telescope HI observations of the ‘Sausage’ cluster to investigate the effect of the 
merger and the shocks on the gas reservoirs fuelling present and future star formation (SF) 
in cluster members. By using spectral stacking, we find that the emission-line galaxies in the 
‘Sausage’ cluster have, on average, as much HI gas as field galaxies (when accounting for the 
fact cluster galaxies are more massive than the field galaxies), contrary to previous studies. 
Since the cluster galaxies are more massive than the field spirals, they may have been able to 
retain their gas during the cluster merger. The large HI reservoirs are expected to be consumed 
within ^0.75 — 1.0 Gyr by the vigorous SF and AGN activity and/or driven out by the out¬ 
flows we observe. We find that the star-formation rate in a large fraction of Ha emission-line 
cluster galaxies correlates well with the radio broad band emission, tracing supernova remnant 
emission. This suggests that the cluster galaxies, all located in post-shock regions, may have 
been underg oing sustained SFR for at least 100 Myr. Th is fully supports the interpretation 
proposed by |Stroe et ah ( |2Q15| ) and Sobral et~ar] ( [2Q15| ) that gas-rich cluster galaxies have 
been triggered to form stars by the passage of the shock. 


Key words: galaxies: active, galaxies: clusters: individual: CIZA J2242.8+5301, shock 
waves, radio continuum: galaxies, radio lines: galaxies 


1 INTRODUCTION 

Galaxy cluster environments have a profound impact on the evo¬ 
lution of cluster galaxies. At low redshifts (z < 0.5) and focusing 
on relaxed clusters, the fraction of galaxies which are star-forming 
drops steeply from field environments, to cluster outskirts and cores 
( |Dressler|1980[|Balogh et al.|1998[[Goto et al.|2003| l . The morpho¬ 
logical transformation of field spirals into cluster ellipticals or SOs 
has been attributed to a number of processes. The dense intracluster 
medium (ICM) could lead to the ram pressure stripping of the gas 
content of field spirals as they accrete onto the cluster (e.g. |Gunn| 
|& Gott|197^|Fumagalli et al.|2014| |. Tidal forces produced by gra¬ 

dients in the cluster gravitational potential or by encounters with 
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other galaxies, can distort infalling galaxies, truncate their halo and 
disk (harassment, [Moore et al.|1996| ) or remove gas contained in 
the galaxy and deposit it into the ICM (strangulation, [Larson et al.j 
[1980j ). All these processes ultimately lead to the removal of gas and 
a truncation of star-formation (SF). 

The effect of relaxed cluster environments on galaxies is evi¬ 
dent using a wide range of diagnostics, which trace different phases 
and time-sc ales of SF. Using UV data produced by young OB stars, 
[Owers et~^ { |2012| ) found galaxies with star-forming trails, which 
they attribute to gas compression by the high-pressure merger envi¬ 
ronment. The UV radiation coming from massive, short-lived stars 
excites emission lines. Lines such as Ha or [OII]3727A probe SF 
on time scales of < 10 Myr. Emission line studies confirm that the 
fraction of star-forming galaxies increases from cluster cores to¬ 
wards field environments (e.g. [Gavazzi et al.[[1998| [Balogh et al.[ 
[T^IPinn et al.||20051 [Sobral et al.||2011| |Darvish et al.||2014t. 

Using far infra-red data (tracing dust obscured SF), [Rawle et al.| 
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^2012^ find that the fraction of dusty star-forming galaxies, com- 
pared to the total star-forming galaxies, increases from low to high 
densities, an effect they attribute to dust stripping and heating pro¬ 
cesses caused by the cluster environment. Similar results are found 
by |Koyama et al.| ( [2013| ). 

Synchrotron emission from supernovae traces SF on longer 
timescales of about 100 Myr ( |Condon|1992| ). Deep radio surveys 
at GHz frequencies indicate that below 100 — 200 jUJy, the number 
of star-forming galaxies dominates over radio-loud active galactic 
nuclei (AGN, e.g. |Padovani et al.|2011|). The num ber of radio-faint 
radio star-forming galaxies ( Morrison & Owen]|2003| ), in clusters 
was found to increase with redshift ( |Morrison|1999| ). Radio broad¬ 
band emission from cluster spirals was also found to be enhanced 
with respect to field counterparts, an effect which can be caused by 
compression of the magnetic fields ( |Gavazzi & Jaffe|1985[|Ander^ 
[sen & Owen|1995| l. 

In addition to probes of past or current SF, CO rotational tran¬ 
sitions can be used as an excellent tracer of molecular gas, which is 
the raw fuel for future SF episodes (e.g. [Leroy et al.|2008| l. Other 
gas phases cannot form stars directly, but they have to cool suffi¬ 
ciently to form cold, dense molecular clouds (see review by jCar-j 
jilli & Walter|2013| ). However, a conversion factor between the CO 
mass and the total molecular gas is needed, which is highly un¬ 
certain (see review by [Bolatto et al.||2013'] l. Instead of using CO 
or other direct tracers of molecular gas, many studies use neutral 
hydrogen HI as a proxy for the molecular gas. Relaxed cluster spi¬ 
rals become increasingly more HI deficient towards cluster cores, 
an effect which does not depend on cluster global properties such 
as X-ray luminosity, temperature, velocity dispersion, richness or 
spiral fracti on (e.g.jMagri et al.|1988{ Cayatte et al.|1990[|Solanes| 


et al.||2008 jCedres et al.||2009[[Hwang & Lee||2009 [Stroe et al.j 


2014a[ Wegner et al.|2015[ [Stroe et al.|2015 ~|Sobral et al.|2015[ l, 


[et al.|200T| [Chung et al.|2009[ l. [0osterloo & van Gorkorn ( [2005[ l 

and Scott et al.[ p012[ l have found galaxies with HI tails, knots 
and filaments, which are possibly caused by ram pressure strip¬ 
ping. Until very recently, HI measurements have been limited to 
the local Universe (z ~ 0). At low redshifts (z ^ 0.06), [Qiengalur[ 
[et al.[ ( |200T] ) studied the A3128 cluster and found that the average 
HI mass for emission-line and late-type cluster members is about 
(8.6 — 8.7) X 10^ Mq. They did not find a statistically significant 
difference between the HI content of emission-line galaxies inside 
and outside the cluster, but the field spirals contain about two times 
more HI than their cluster counterparts. Pioneering work by [Verhei^ 
[jen et al.[ ( [2007[ ), [Lah et al.[ ( [2007[ ) and [Lah et al.[ ( [2009] ) used direct 
detections and stacking to measure the HI content of cluster galax¬ 
ies up to z ~ 0.4. [Verheijen et al.] ( [2007[ ) surveyed two z ~ 0.2 clus¬ 
ters with very different morphologies: the relaxed, massive galaxy 
cluster A963 and the low-mass, diffuse cluster A2192. They detect 
only one HI galaxy within 1 Mpc of the centre of each cluster. By 
stacking galaxies with known redshifts, they make a clear detection 
of HI for blue galaxies outside the clusters, but no such detection 
was made for cluster galaxies. In a detailed study of the cluster 
A370 at z ^ 0.37, [Lah et al.[ ( |2007] ) use spectral stacking to mea¬ 
sure HI in gas-rich galaxies lying outside or at the outskirts of the 
cluster. 

As discussed previously, relaxed cluster environments are be¬ 
lieved to suppress SF by removing cold gas form their host galax¬ 
ies. At z < 0.3 between 10 — 20 per cent of clusters are undergoing 
mergers ( [Katayama et al.|2003l [Sanderson et al.|2009[ [Hudson et[ 

[al.|2010[l and this fraction is ex pected to increase steeply beyond 

z = 0.4 ( [Mann & Ebeling|2012[ ). The effect of cluster mergers on 
the SF activity and gas content of galaxies is disputed. Most studies 
find that cluster mergers trigger SF ( Miller & Owen|20(}^ [Umeda 


but a few studies find they quench it (e.g. Poggianti et al. 12001) or 

that they have no direct effect (e.g. [Chung et al.[2010[ l. 

An interesting subset of clusters are those hosting radio relics, 
extended patches of diffuse radio emission tracing merger-induced 
shocks (e.g. [Ensslin et al.|1998[l. The Ha properties of radio-relic 
clusters Abell 521, CIZA 12242.8-^5301 and IRXS 10603.3-^4214 
( [Umeda et al.[[2004l [Stroe et ^[2014a[ [2015[ ) indicate that the 
merger and the passage of the shocks lead to a steep SF increase 
for < 0.5 Gyr. The fast consumption of the gas ultimately quenches 
the galaxies within a few hundred Myr timescales ([Roediger et al.[ 
[20T4l ). 

The Ha studies of [Umeda et al.[ ( [2004] ), [Stroe et al.[ ( [2014a[ ) 
and [Stroe et al.[ ( |20T5] ) are tracing instantaneous (averaged over 10 
Myr) SF and little is known about SF on longer timescales and 
the reservoir of gas that would enable future SF. An excellent test 
case for studying the gas content of galaxies within merging clus¬ 
ters with shocks is CIZA J2242.8-1-5301 ( [Kocevski et al.[[2007] l. 
For this particular cluster unfortunately, its location in the Galac¬ 
tic plane, prohibits studies of the rest-frame UV or FIR tracing 
SF on longer timescales, as the emission is dominated by Milky 
Way dust. However, the rich multi-wavelength data available for 
the cluster give us an unprecedented detailed view on the inter¬ 
action of their shock systems with the member galaxies. CIZA 
J2242.8-1-53 01 is an extremely massive ( [Jee et al.][2015[ [^awson 


et al.|2015" M 200 ^ 2 x lO^^M©) and X-ray disturbed cluster \A~ 


matsu & Kawahara[[2013| [Ogrean et al.[[2013[ [2014[ ) which most 


likely resulted from a head-on collision of two, equal-mass systems 
( [van Weeren et al.|2011[ [Dawson et al.|2015[ l. The cluster merger 
induced relatively strong shocks, which travelled through the ICM, 
accelerated particles to produce relics towards the north and south 
of the cluster ( [van Weeren et al.|2010[ [Stroe et al.|20T3] ). There is 
evidence for a few additional smaller shock fronts throughout the 
cluster volume ( [Stroe et al. [20 13 [[Ogrean et al.[2014[ ). Of particular 
interest is the northern relic, which earned the cluster the nickname 
‘Sausage’. The relic, tracing a shock of Mach number M ^3 ( [Stroe[ 
[etal.[2014c] l, is detected over a spatial extent of ^ 1.5 Mpc in length 
and up to ^ 150 kpc in width and over a wide radio frequency range 
(150 MHz —16 GHz; [Stroe et al.|20"T3l[2014b] l. There is evidence 
that the merger and the shocks shape the evolution of cluster galax¬ 
ies. The radio jets are bent into a head-tail morphology aligned with 
the merger axis of the cluster. This is probably ram pressure caused 
by the relative motion of galaxies with respect to the ICM ( [Stroe[ 
[et akpOl^j l. The cluster was also found to host a high fraction of 
Ha emitting galaxies ( [Stroe et al.|2014a[|2015[ ). The cluster galax¬ 
ies not only exhibit increased SF and AGN activity compared to 
their field counterparts, but are also more massive, more metal rich 
and show evidence for outfiows likely driven by super-novae (SN) 
( [Sobral et al.[[2015[ l. [Stroe et al.[ ( [2()15[ ) and [Sobral et al.[ ( |2015[ ) 
suggest that these relative massive galaxies (stellar masses of up 
to ~ Mq) retained the metal-rich gas, which was trig¬ 

gered to collapse into dense star-forming clouds by th e passage of 
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the shocks, travelling at speeds up to ~ 2500 km s“^ ( [Stroe et ah 
[2014c[ l, in line with simulations by [Roediger et al.[ ( [2014[ ). 

In this paper we focus on the effect of the massive cluster 
merger and travelling shocks in the ‘Sausage’ cluster on the HI 
content of the galaxies, tracing the gas that may fuel future SF. 
We place this in context of other phases of SF, averaged over 
short timescales (^10 Myr, Ha data) and averaged over longer 
timescales ('^ 100 Myr, radio broad band data) SF episodes. 

The structure of the paper is as follows: in ^we discuss the 
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Table 1. Details of the JVLA observations of the ‘Sausage’ cluster taken in L band (1.5 GHz), combining all four configurations of the telescope for a total of 
> 26 h of observing time. 



L-band A-array 

L-band B-array 

L-band C-array 

L-band D-array 

Observation dates 

Total used on source time (h) 
Integration time (s) 

May 11,2014 
~6.5 

1 

Oct 31, 2013 
-6.5 

3 

Sep 2, 2013; Jul 3, 2013 
- 10 

5 

Feb 3, 2013; Jan 31, 2013; Jan 27, 2013 
-3.5 

5 


observations and the reduction of the HI, optical and broad band 
data; in ^we discuss the HI detections, stacking, masses and how 
these correlate with Ha and radio luminosities; in ^we discuss the 
implication of our results for future SF episodes and compare them 
with other HI studies of clusters. Finally, ^presents a summary of 
the results, placing them in context of the SF history of the cluster. 
At the redshift of the ‘Sausage’ cluster, z ^ 0.188,1 arcsec covers a 
physical scale of 3.18 kpc and the luminosity distance is Jp ~ 940 
Mpc. All coordinates are in the J2000 coordinate system. We use a 
|Chabrier| ( |2003| ) initial mass function (IMF) throughout the paper. 
We correct measurements from other papers accordingly. 


2 OBSERVATIONS & DATA REDUCTION 

In our analysis, we combine radio spectral line data tracing HI, 
broad-band radio data, optical imaging and spectroscopy of passive 
and star-forming galaxies in and around the cluster. 

2.1 HI data 

The ‘Sausage’ cluster was observed with the Westerbork Synthesis 
Radio Telescop^in the maxi-short configuration during the second 
half of 2012, for a total of 26 12-h tracks. The observations were 
taken in 4, slightly-overlapping bands of 10 MHz each and central 
frequencies 1180, 1188, 1995 and 1203 MHz, respectively, there¬ 
fore fully covering the 1175 — 1208 MHz range. Re-circulation was 
used to have 512 channels per band (with XX and YY polarisations 
only). Sources CTD93 and 3C 147 were used as flux calibrators. 

The velocity resolution (after Hanning smoothing), at the red- 
shift of the cluster, is about 20 km s“^, sampled a channel width of 
^ 9.9 km s“^ The velocity range covered is about 4700 km s“^, 
corresponding to an HI redshift range of 0.184 — 0.199. The red- 
shift range covers the cluster volume within — Icr to +3(7 of the 
cluster redshift z = 0.188, where a = 0.04 is the cluster velocity 
dispersion ( [Dawson et al.|2015] l. The HI observations cover well 
the distribution of the Ha emitting galaxies (z = 0.190 + 0.010; 
jSobral et al.|2015| ). 

Significant radio frequency interference (RFI) is known to be 
present at the WSRT at frequencies covering the HI redshift range 
of 0.1 to 0.25, caused by geo-positional systems such as GPS and 
GLONASS. However, at the time of our observations, the frequen¬ 
cies corresponding to the redshift of the cluster (1185 — 1200 MHz) 
were still fairly free of such RFI. With the foreseen deployment 
of the GALILEO geo-positional system, the situation in this fre¬ 
quency range will worsen dramatically in the near future. 

To remove any residual RFI, we performed our RFI flagging 
on the Stokes Q (i.e. XX-YY) component on the data. Given the po¬ 
larised nature of RFI, this removed most of the astronomical signal, 

^ http://WWW.astron.nl/radio-observatory/astronomers/ 
observing-wsrt/observing-wsrt 


but left the RFI mostly intact. The flags found for Stokes Q were 
then applied to original XX and YY visibilities. Moreover, because 
the RFI in our data is broad band, before flagging we performed a 
smoothing in frequency to enhance the ‘sensitivity’ for RFI. These 
procedures worked very well and the final data cubes do not show 
any effects of residual RFI, while the noise level (75 /iJy beam“^ 
over 20 km s“^) is very close to what is expected for the integration 
time. 

Once the data were flagged, we applied standard calibration 
procedures to the data using the Miriad software ( jSault et al.|1995j ). 
The continuum emission was removed from the data by fitting, to 
all channels, a 3^^ order polynomial to each visibility spectrum 
(‘uvlin’). The choice of the order depends on out to which ra¬ 
dius there are significant continuum sources and on bandwidth. The 
higher the order, the better the sources are removed, however the 
noise, after subtraction, increases. A 3^^ order polynomial fit rep¬ 
resents a compromise between sufficient removal of the continuum 
and little increase in the noise level of the line cube, as shown by 
jSaultj l p^^ . Because this does not take into account the presence 
of any possible HI emission, the spectra of individual detections, 
and of the stacked spectra we discuss below, are corrected for this 
over-subtraction (see Section [33] ). 

The synthesised beam of the WSRT observations is 24.9 x 
18.5 arcsec^ at a position angle of 165.8°, or 79.3 x 58.9 kpc^. 

2.2 Optical imaging and spectroscopy 

In order to locate the spatial and velocity position of the HI signal, 
we use the Keck and William Herschel Telescope spectroscopy. 
Data of both passive and star-forming galaxies in the field of the 
cluster were presented in [Sobral et al.[ ( [20l3] ), [Dawson et al.[ ( [2015[ ) 
an d [Stroe et'^ ( [2015[ ). 

Galaxies are categorised as passive or emission-line based on 
spectral features. Emission line galaxies were selected through the 
presence of the Ha emission-line (with a Ha-i-[NII] rest-frame 
equivalent width larger than 13A), tracing hot ionised gas, indicat¬ 
ing the presence of SE and/or radio-quiet AGN (broad and narrow 
line AGN). The spectra of passive galaxies display Baimer absorp¬ 
tion features and/or no Balmer series emission lines. The passive 
galaxies are undetected in Ha at the 13 A level equivalent width 
level. The sample was divided into three categories: 

(i) passive galaxies inside the cluster, 

(ii) emission-line galaxies within the cluster, 

(iii) emission-line galaxies in the field around the cluster. 

The number of sources in each sub-sample can be found in Tablej^ 
The emission-line galaxy sample is dominated by star-forming 
galaxies with a ^ 20 — 30 per cent contribution from broad and 
narrow line optical AGN (see Table [Sobral et al.[[2015[ ). The 
cluster members were chosen to be located at a projected distance 
of less than 1.85 Mpc away from the cluster ‘centre’, in line with 
the definition of cluster membership from [Stroe et al.[ ( |T015[ ). Line- 
emission galaxies outside the cluster were defined to lie outside 
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Table 2. Number of sources with spectroscopic data, separated by galaxy type. Not all spectroscopic sources are covered by the redshift range of the HI 
observations. 


Sample Total number Sources within HI z range Reference 


Emission line, field 

39 

22 

Sobral et al. 

2015 

';S 

troe et al. 

2015 

Emission line, cluster 

54 

45 

Sobral et al. 

2015 


troe et al. 

2015 

Cluster star-forming 

48 

39 


Sol 

3ral et 

al. 




Cluster AGN 

6 

6 


Sobral et al. 

i2015 



All emission-line 

93 

67 

Sobral et 

W:(|2015|;S 

troe et 

al. 

j2015 

Passive, cluster 

184 

154 

Dawson et al 



All 277 221 


of the 1.85 Mpc radius. The redshift distribution of the galaxies is 
plotted in FigureWe note that, as shown by |Sobral et al.| ( [20T5] ), 
the samples of cluster and field line emitters are selected uniformly, 
down to a similar star-formation rate (SFR). 

The spectroscopy is supplemented with Subaru, Canada- 
France-Hawaii, William Herschel Telescope and Isaac Newton 
Telescope broad band (BB) and narrow-band (NB) imaging trac¬ 
ing the Ha emission-line in galaxies at the cluster redshift ( |Stroe| 
|et al.|2015] l. Ha luminosities for each source were calculated using 
the excess of NB emission compared to the i BB filter (for method 
and details see |Stroe et al.|2015] l. In the analysis, we also employ 
stellar masses derived using the method described in|Sobral et al.| 
f20T5\ . 


2.3 Broad-band radio data 

We identified radio counterparts to the optical sources by cross¬ 
matching with a deep, high resolution 1.5 arcsec) image of the 
cluster, centred at 1.5 GHz, produced using the upgraded Jansky 
Very Large Array. 

Deep JVLA observations of the cluster were taken in the 1 — 2 
GHz L-band in A, B, C, and D-array configurations. An overview 
of the observations is given the Table In total, 16 spectral win¬ 
dows with 64 channels, each covering 64 MHz of bandwidth, were 
recorded. The data reduction for each observing run was carried out 
separately, using CAS^version 4.2. 

As a first step, the data was Hanning-smoothed and correc¬ 
tions for the antenna positions and elevation dependent gains were 
applied. We then obtained an approximate bandpass solution us¬ 
ing observations on the primary calibrators (3C147, 3C138). We 
applied the bandpass solutions to the data and flagged RFI in an 
automated way using the AOFlagger ( [Offringa et al.|2010] ). This 
initial bandpass correction was performed to avoid flagging of data 
due to the bandpass variations across the spectral windows. Af¬ 
ter flagging, we obtained gain corrections on the primary calibra¬ 
tors using 5 channels centred at channel 30. These gain solutions 
were obtained to remove the time-varying gains. We pre-applied 
these solutions to find delay terms and bandpass solutions. We then 
re-determined the gain solutions but now using the full channel 
range pre-applying the bandpass and delay corrections. We then 
pre-applied these solutions to obtain gain solutions on our sec¬ 
ondary calibrator J0542-I-4951. The cross-hand delays were solved 
for using the calibrator 3C 138. The channel dependent polarization 
leakage terms and polarization angles were set using 3C 147 and 
3C 138, respectively. For observing runs longer than four hours the 


^ http://casa.nrao.edu 


leakage terms were determined from scans on the secondary cali- 
bratorJ0542+495]|!] 

We bootstrapped the flux-scale from our primary calibrator 
observations to find the flux-density of J0542-I-4951. The flux-scale 
was set using the default settings of the task set jy. As a final step, 
the calibration tables were applied to the target data. We averaged 
the target field data by a factor of 3 in time and frequency, to reduce 
the data volume for imaging. 

In the next steps, the calibration solutions were refined us¬ 
ing several cycles of phase-only and amplitude and phase self¬ 
calibration. For the imaging we employed w-projection (|Cornwell| 
|enil|2008l [2005} and MS-MFS clean ( |Rau & Comwell||2011t 
with nterms=3. Clean boxes were set by running the Python Blob 
Detection and Source Measurement (PyBDSlQ). A few additional 
clean regions were added manually for diffuse sources. 

After self-calibrating the individual datasets, we combined all 
the datasets to make one deep image. Two more rounds of ampli¬ 
tude and phase self-calibration (on a 1 h timescale) were carried 
out on the combined dataset. During the self-calibration the am¬ 
plitude scale was allowed to drift freely. This was needed to fully 
align the different datasets and spectral windows and avoid strong 
artefacts around a few bright sources located in the field of view 
(FOv|^ We made a final image using Briggs (robust=0) weight¬ 
ing and corrected for the primary beam attenuation. We checked the 
flux-scale of the image against our previous 1.38 GHz WSRT ob¬ 
servations of the cluster ( [van Weeren et al.|2010[[Stroe et al.|20T3] l. 
This was done by checking the integrated fluxes of 10 sources and 
scaling the fluxes from 1.4 to 1.5 GHz assuming a spectral index 
of —0.7, the canonical spectral index for bright radio sources (e.g. 
|Condon|1992| ). The JVLA fluxes were divided by a correction fac¬ 
tor of 1.4 ±0.1 to re-align the flux-scale to the WSRT scales. The 
resolution of the VLA image is 1.5 arcsec x 1.4 arcsec, with a po¬ 
sition angle of 86.5 degrees. 

3 RESULTS 

3.1 Direct detections 

Six galaxies are tentatively detected in HI. The detection criterion is 
signal above 5(7 over at least two consecutive channels. Tablej^lists 
the redshifts and HI masses of the direct detections. The narrow HI 
profiles (^ 40 km s“^) of the directly detected sources indicates, if 

^ The polarization results will be discussed in a forthcoming paper (van 
Weeren et al. in prep.) 

^ http://dl.dropboxusercontent.eom/u/1948170/html/index.html 
^ In principle, the global amplitude scale could have been preserved, but 
such an option is not offered in CASA at the moment. 
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Figure 1. Redshift distribution of sources in each galaxy sample (emission¬ 
line cluster members, emission-line field galaxies and cluster passive galax¬ 
ies). The HI redshift coverage of the WSRT data is marked by the vertical 
dashed lines. 

Table 3. Properties of the direct detections. The positional uncertainties are 
15 arcsec, the error in the redshifts is 0.0005 and the error in the HI mass 
about 0.4 • 10^ M©. The stellar masses reported are for the closest face-on 
spiral found within 7 arcsec. 


RA 
(h m s) 

DEC 

(O , 

z 

Mm 

(109 Mq) 

M* 

(109 Mq) 

Mhi/M* 

22 41 51.1 

52 52 55 

0.18893 

2.1 

2.0 

1.0 

22 42 56.5 

52 57 21 

0.18916 

2.6 


- 

22 43 14.3 

53 04 57 

0.18536 

1.2 

9.7 

0.1 

22 43 23.2 

53 04 41 

0.18536 

2.2 

** 

- 

22 41 30.4 

53 05 58 

0.18486 

1.7 


- 

22 43 43.7 

53 09 43 

0.18994 

1.6 

7.9 

0.2 


* A bright star overlaps the position of the galaxy, so no reliable 
counterpart can be found. 

* No face-on spiral was found nearby the HI detection indicating the 
detection is spurious. 

real, they are most probably oriented in the plane of the sky. For a 
given HI mass, the peak flux for a face on galaxy is higher than for 
an edge-on galaxy because the same amount of flux divided up in 
fewer channels. Therefore, we are biased towards detecting face-on 
galaxies. However, given the very narrow profiles of these tentative 
detections, they could also be high noise peaks. 

We note that none of the six HI direct detections have counter¬ 
parts (matches within 10 arcsec) in our spectroscopic data or in the 
Ha catalogue. This indicates that the sources have faint Ha emis¬ 
sion, below the equivalent width detection threshold (13A). Their 
SFR are therefore below ^0.35 Mq yr“^. 

Given the positional accuracy coupled with the large beam of 
WSRT finding the right optical counterpart is very challenging. A 
few (1—6) potential optical hosts are found for the tentative HI di¬ 
rect detections, but most sources are unresolved and faint (/ band 
magnitude on average fainter than 20). Therefore reliable photo¬ 
metric redshifts could not be derived and we cannot confirm them 
as being located at z ~ 0.2. However, even if these sources were 
z ^ 0.2 galaxies, they would have small stellar masses. For exam¬ 
ple, using the closest optical match as galaxy host, we And that 
their stellar masses are very small (< 1 x lO^M©). This would im¬ 
ply unrealistically high gas to stellar mass ratios. Additionally, the 
morphology of these close optical sources does not match face- 
on spiral galaxies, as we expect. These arguments indicate that the 
closest sources are not the correct matches. For three out of the six 
HI detections there are large, face-on spiral galaxies in their vicin¬ 
ity (within 7 arcsec), which could be the correct optical counter¬ 
part. These three galaxies have stellar masses of ^ (2 — 10) x 10^ 
Mq, which indicates 0.1 — 1 atomic gas to stellar fractions. For two 
detections only small sources are located in the vicinity and no ob¬ 
vious face-on galaxies are found nearby the HI, indicating these are 


spurious detections (noise peaks in adjacent channels mimicking a 
signal). In one case, a bright star located at the location of the HI 
detection prevent correct optical identification. 


3.2 HI stacking 


Since only six galaxies were directly detected, we use spectral 
stacking to measure the average HI content of the galaxy samples. 

We use the optical positions to extract radio spectra for each 
galaxy, summing the flux within an elliptical aperture equal to the 
FWHM of the synthesised beam (24.9 x 18.5 arcsec^). This corre¬ 
sponds to a spatial scale of 79.3 x 58.8 kpc^, matched to the phys¬ 
ical size of the galaxies, which at z ~ 0.19 are unresolved in the 
HI observations ( [Verheijen et al.|2007]l (also in line with the galaxy 
sizes presented in |Stroe et al.|2014a 120151 1. We further tested the 
effect of the aperture size on the final stack using apertures of sizes 
ranging from 0.1 of the FWHM up to 2 times the FWHM. 

To test the effect of aperture size on the final HI stack, we use 
elliptical apertures in size equal to a fraction of the FWHM on both 
the width and height of the synthesised beam. We use apertures 
from 0.1 to 2 times the FWHM of the radio beam. We follow the 
procedure described in § |3.2| to extract spectra at the positions of 
the cluster emission-line members. We And that the peak of the HI 
detection remains relatively stable if the aperture is at least 0.9 of 
the FWHM size (see top panel. Figure]^. 

For each aperture size (from 0.1 to 2.0 times the FWHM of 
the radio beam, see Figure [^, we also measure the HI mass in the 
way described in §3.4| As shown in the bottom panel of Figure 
we And that the HI mass is relatively stable as function of aperture, 
but it peaks when the size of the aperture is equal to the FWHM 
of the synthesised beam. The FWHM size of WSRT is also well 
matched to the expected HI disk size of galaxies at z ~ 0.19 ( |Ver-| 
[heijen et al.|2007| l. We therefore choose to measure HI quantities 
within apertures equal to the FWHM of the synthesised beam. 

We extract spectra for all galaxies whose redshift falls within 
the HI redshift range probed by our WSRT observations (0.184 < 
z < 0.199). To study the noise properties, we extract spectra in 
sky positions shifted by —60 arcsec in the RA direction, but us¬ 
ing the same redshifts as the sample of galaxies. This method cap¬ 
tures the effect of increasing noise towards the edges of the band¬ 
pass. Note however that the noise could be overestimated. Given the 
large source density of the cluster held, a shift in sky position does 
not guarantee we will be measuring pure noise, but some apertures 
could partially fall on undetected source. 

Before stacking, we correct the galaxy and noise spectra for 
the WSRT primary beam, which is a function of distance from the 
pointing centre and observing frequency: 

A(r, v) = cos^(cvr) (1) 


where c = 68 is a constant, r is the distance from the pointing center 
in degrees and v is the observing frequency in GHz. By correcting 
for the primary beam, we account for the effect of noise increasing 
towards the FOV edges in both the galaxy and their associated noise 
spectra. 

The spectra are then shifted to the rest-frame velocity using 
their spectroscopic redshift. We use the radio deflnition of veloc¬ 
ity: F = c(vhi - Vobserved)/vHT where c is the speed of light, Vri 
is the restframe frequency of HI and Vobserved is the observing fre¬ 
quency. The galaxy spectra and the noise spectra are co-added us¬ 
ing a weight based on the primary beam (oc A^(r, v)), which ac¬ 
counts for the fact that the noise for galaxies away from the held 
centre is larger in proportion to the primary beam. To obtain the 
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Figure 2. Top: HI peak emission flux in the cluster line emitter stack as 
function of the aperture size used for extracting the individual spectra. 
The horizontal scale indicates the fraction by which the aperture width and 
height are scaled as a function of the FWHM. The trend in the peak emis¬ 
sion does not change before and after correcting for the continuum over¬ 
subtraction (see §3.3} . Bottom: Average HI mass of cluster line emitters 
(see §3.4} , as function of the aperture size used for extracting the individual 
spectra. The horizontal scale indicates the fraction by which the aperture 
width and height are scaled as a function of the FWHM. 

correct flux density scale, we normalise the stacked spectrum by 
the integral of the synthesised beam, integrated over the same spa¬ 
tial region used for extracting the galaxy spectra. 

After stacking and normalising the spectra, we filter the data 
using a second-order Savitzky-Golay filter (SG; [Savitzky & Go-| 
|lay|1964} , which convolves the data with a polynomial filter. Given 
that our line profiles are resolved (FWHM of the stacks with detec¬ 
tion are >150 km s“^, compared to a channel width of > 20 km 
s“^ see Table 1^, the method reduces the noise, while preserving 
line profiles (see for example [Morabito et al.|2014] ). We tested the 
method using different window sizes. The effect of the SG filtering 
with increasing window lengths is shown in Figure Out of the 
window size tested, we finally choose a filter window of ~ 300 km 
s“^, which provides minimal noise, while preserving the width and 
height of the signal, thus maximising the signal to noise of the pos¬ 
sible HI detection peaks. We also tested other smoothing kernels 
(e.g. moving boxcar), which resulted is similar results, but with a 
widening of the profile and reduction of the peak strength. 

Separate stacks are produced for the sample of passive cluster 
galaxies, line emission cluster members and line emission galaxies 
located within the field environment around the cluster. Line emit¬ 
ters include both star-forming galaxies and optical AGN (see also 
§ |2.2} . Finally, we produce a master stack of all the galaxies avail¬ 
able. The number of galaxies per each velocity channel for each 
stack is shown in Figure and Table The number of sources 
in the stack naturally peaks at the 0-velocity position, but dwin¬ 
dles towards higher relative velocities. This effect is governed by 
where the redshift of each source falls within the WSRT HI band¬ 
pass. Due to extensive spectroscopy from Keck aimed at obtaining 
a dynamical analysis of the cluster that specifically targeted the red 
sequence, the number of passive cluster galaxies far outnumbers 
the number of emission-line galaxies. 

The asymmetry (about the central position) in the number of 



Figure 3. Noise reducing properties of the Savitzky-Golay filter with in¬ 
creasing filter window width demonstrated on the cluster member galaxy 
stack. Colours of the lines go from lighter to darker shades in increasing 
window length. The black line at the foreground of the figure is smoother 
spectrum after applying the final window length choice, ~ 300 km s“^ 
which is also the one we finally applied to the data. Similar results are ob¬ 
tained for the galaxy stacks, as well as their associated noise stacks. 


sources for which data exists, especially visible in the passive clus¬ 
ter galaxy sample, is caused by the discrepancy in the nominal red- 
shift of the cluster (z = 0.188, recently derived from more than 
200 spectra by Dawson et al. 1120151 1 and the outdated z = 0.192 
( |Kocevski et al.||2007} which was used for creating the WSRT 
HI setup. Our spectroscopic measurements are therefore biased to¬ 
wards lower redshifts, compared to the HI coverage of the WSRT 
data (see §2.2| and Figure [^. If a redshift of a source falls at the 
middle of the WSRT HI band coverage, the frequency coverage 
is symmetric about the observed frequency of the HI. However, a 
lower redshift (than the central redshift) is equivalent to a source 
having a wider frequency coverage at frequencies lower than the 
HI and a narrower coverage at higher frequencies. When translat¬ 
ing to a rest-frame velocity, there is a preferential data coverage of 
the positive restframe velocities. The missing data at larger abso¬ 
lute restframe velocities drives the noise to higher values in those 
regions. 


3.3 Correcting for the over-subtracted continuum emission 

As mentioned previously, the bulk of our galaxies do not have a 
direct detection of the HI line. Hence, the continuum emission sub¬ 
traction (see could not take into account the presence of HI 
emission and leads to an over-subtraction of the continuum where 
the putative HI line is located. The over-subtraction is not visi¬ 
ble (and relevant) in the individual spectra, but it shows up in the 
stacked spectrum. As expected, the over-correction of the contin¬ 
uum is evident in the case of the emission-line galaxy samples, 
where the HI signal is located on top of a broad, negative dip. No 
evident negative trough is present in the passive member sample, 
where less HI is expected. 

We apply a two step process to correct for the over-subtraction 
of the continuum. Firstly, to measure the possible extent of the HI 
in the Ha galaxies, we fit Gaussian profiles to their stacked profiles. 
We select data at least dcJoauss away from the peak of the Gaussian, 
where CJoauss is the dispersion of Gaussian profile. Additionally, we 
discard the data at relative restframe velocities higher than 1000 
km s“l. These cuts are employed to exclude any HI signal from 
the estimation of the continuum, but not include very noise edge 
channels, where only a few galaxies are stacked (see Figure]^. In 
the case of the passive members, a similar procedure is applied, 
but since we do not have a clear detection of emission, the velocity 
range between 250 and 1000 km s“^ is used. 

We fit the channels free of HI with a 2nd order polynomial and 
subtract the fit from the data to correct for the over-subtracted con- 
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Emission line, field All emission line All 

Emission line, cluster Passive, cluster 



Figure 4. Number of sources stacked in each velocity channel, for the sep¬ 
arate stacks. The passive cluster galaxy stack (purple line) benefits from a 
factor of > 5 more sources than the stacks for the Ha line emitters (mostly 
star-forming and some AGN) galaxies inside (green line) and outside the 
cluster (black line). Note how the number of sources dwindles towards the 
edges of the WSRT HI bandpass. The asymmetry in the number of sources 
is caused by the redshift distribution of the stacked sources (see Figure [T], 
preferentially located at redshifts lower than the central WSRT HI coverage. 

tinuum (see example in Figure [^. The correction is substantial for 
the cluster Ha line galaxies, but the results indicate that the contin¬ 
uum has not been significantly over-subtracted for passive cluster 
members and those line emitters located in the field environment. 


3.4 Measuring the HI signal and its significance 

After the stacked spectra are filtered and corrected for the con¬ 
tinuum over-subtraction, we fit Gaussian profiles around the 0- 
velocity position. Results are shown in Figure The parameters 
of the Gaussian fits can be found in Table|4] We calculate the RMS 
from the noise stack. Note that even though the number of galaxies 
in the cluster and field stacks is similar, the noise levels achieved are 
a factor of ~ 4 different. This is because field galaxies are preferen¬ 
tially located at large distances away from the FOV centre, meaning 
that the noise levels at their location are higher. 

In the case of the Ha emission-line galaxies, a clear peak is 
found around the 0-velocity position and a Gaussian profile could 
clearly be fit. In the cluster line emitter stack, we reach a noise level 
Orms of 5.9 jUJy and detect HI at a peak significance of 7.2(7. The 
results are virtually unchanged if the 6 AGN are removed from the 
stack. No clear detection of HI is made in the case of the passive 
galaxies, despite the much larger number (2.5 times) of galaxies 
stacked, as compared to the line emitter sample (HI measured with 
a peak of 5.3 jUJy, with a cjrms = 2.1 jUJy). The velocity posi¬ 
tion of the putative peak (~ —100 km s“^), highly offset from 0, 
is likely a spurious peak and also indicative of a non-detection of 
HI. The high offset is highly unlikely to be caused by stripping, 
given we are averaging across 150 passive galaxies within the clus¬ 
ter, which have random motions in the cluster potential. In some of 
the stacks (e.g. emission cluster galaxy stack and the emission line 
field galaxy stack), there are additional peaks off-centred from the 
0-velocity. In theory, additional peaks in cluster stack, for example, 
could be caused by tidally stripped tails, pointing in the same red- 
shift direction, such that they add coherently when stacked. How¬ 
ever, this is highly unlikely as the cluster galaxies are moving on a 
range of orbits within the cluster potential. Even if HI tails existed, 
they would have a variety of orientations. When stacking the galax¬ 
ies, the tails would therefore not add coherently. Therefore, we be¬ 
lieve these to peaks to be caused by noise variations and low-level 
systematics. Note that at higher restframe velocities, the number 
of sources for which data exists dwindles (as shown in Figure]^. 


For example, at ±600 km s“^, the number of sources with data in 
that velocity channel is already half that at velocity 0. That means 
that the noise at larger restframe velocities will be higher than the 
noise around the 0-velocity. Therefore the significance of off-center 
peaks is actually very low (at least a factor \/2 lower than if located 
at 0-velocity). 

The stack using all the galaxies reaches a very low noise level 
of 3.1 jUJy. The HI signal peaks at a significance of 4.6(7 rms^ value 
mainly due to the contribution of SF and AGN galaxies. 

3.5 HI masses 

We use the following relation to convert from radio flux S into HI 
mass (MHi; |Wentzel & van Woerden|1959[|Roberts|196^ : 


Mhi 

236 


f fSvdV \ 

Mq 

TT^ ' 

i^Mpcy ' 

^mJy km s“^ J 


where Mq is the mass of the Sun, z = 0.192 is the mean redshift of 
the sample of galaxies, Dl = 940.7 Mpc is the luminosity distance 
at that redshift and / SydV is the average of the HI emission over 
a restframe velocity range. As mentioned in §3.2| the velocity is 
defined as V = (^(Vhi — l^observed)/^HT 

For the stacks with HI detections (Figure [^, the HI mass is 
averaged over 2(7Gauss (2 times the Gaussian dispersion) on either 
side of the peak position. For the passive population, we use the 
range within 200 km s“^ from 0-velocity position. The error in 
the HI mass is calculated by propagating the RMS error through 
equation 1^ using the same velocity range used for the integration 
of the signal. 

Note that even though we use the filtered stacks to calculate 
the HI masses, similar results would be obtained even if the orig¬ 
inal data is used. Although the quality of the spectra improve, the 
errors on the final HI mass do not heavily depend on the filtering. 
As shown in equation]^ the mass is effectively an integral over the 
profile. Hence, the effects of SG filtering are reduced because av¬ 
eraging over a velocity range equivalent to smoothing down to a 
resolution equal to that velocity range. 

We find that the average HI mass for the emission-line clus¬ 
ter galaxies (Mhi = (2.50 ± 0.62) x 10^ Mq) is a factor of 1.3 
higher than the mass of cold neutral gas in their field counterparts 
(Mhi = (1.86 ± 1.20) x 10^ Mq). However the difference is not sig¬ 
nificant. Additionally, the cluster galaxies are on average ^1.5 
times more massive than their field counterparts (see also Table [^. 
The average stellar masses (M^) are calculated over the same galax¬ 
ies stacked for the HI analysis and the error reported is the standard 
deviation of the sample. Therefore, when accounting for the differ¬ 
ences in stellar masses between the cluster and the field galaxies, 
the cluster line emitters are consistent with being as gas rich as the 
field counterparts. The fraction of neutral atomic gas to stellar mass 
Mhi/M^ is 0.39 ± 0.26 and 0.34 ± 0.09 for field and cluster galax¬ 
ies, respectively. This is a surprising result, as cluster galaxies have 
been found to contain less HI than galaxies in the field (e.g.|Solanes| 
let aL|2001| ). 

The HI masses of the directly detected sources (Table are 
in the same range as the values for the average emission-line stacks 
(Table 1^. However, we do not directly detect any of the emission¬ 
line systems. This indicates that we are biased towards detecting 
face-on sources, while the galaxies used for stacking have random 
orientation (with the bulk being oblique to edge on). 

Within our noise limits, we do not detect any significant 
HI within the passive cluster galaxy population (Mhi = (0.21 ± 
0.15) X \0^Mq). The passive population has at least nine times less 
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Table 4. Peak number of sources in every stack created for the different galaxy samples, the RMS noise value obtained for each stack and the parameters of 
the Gaussian fit to the HI signal. HI emission is securely detected, at its peak, in the emission-line galaxy and all galaxy stacks, but not in the cluster passive 
galaxies. When integrating the HI signal to estimate an HI mass, we obtain clear detections for the line-emitters. 


Sample 

Number 

C^RMS 

(MJy) 

HI peak 

(MJy) 

Peak significance 
Orms 

HI velocity 
(km s“^) 

HI width 
(km s“^) 

Mm 

(10® Mq) 

(109 Mq) 

Mhi/M* 

Emission line, field 

22 

21.3 

60.4 

2.8 

1.0 

69.3 

1.86 ±1.20 

4.8±0.8 

0.39 ±0.26 

Emission line, cluster 

45 

5.9 

42.7 

7.2 

4.6 

137.0 

2.50 ±0.62 

7.4±0.5 

0.34 ±0.09 

All emission-line 

67 

8.2 

45.1 

5.5 

-0.7 

104.7 

2.00 ±0.67 



Passive, cluster 

154 

2.1 

5.3 

2.5 

-97.4 

90.1 

0.21 ±0.15 

25.6 ±0.4 

<0.008 ±0.006 

All 

221 

3.1 

14.3 

4.6 

-29.5 

107.0 

0.60 ±0.30 




fin 


Noise stack 

Galaxy stack 

_ 

Gaussian fit 







(a) Emission line cluster galaxies. 



(b) Emission line field galaxies. 


(c) All emission-line galaxies. 




(d) Cluster passive galaxies. 


(e) All galaxies. 


Figure 5. Einal HI stacks after filtering and correcting for continuum over-subtraction. The solid black lines show the galaxy stacks and the gray lines show 
the corresponding noise properties. The thin, solid red lines show the Gaussian fits to the profiles located around a 0 restframe velocity. The vertical, dashed 
red line show the position of the peak, which is labelled with its significance based on the underlying noise properties. The RMS noise is calculated using only 
the data between the vertical, dashed gray lines, where at least 80% of the peak number of sources are stacked. Note clear detections of HI are made for the 
line emitters, while no detections are made in the case of the cluster passive galaxies. 
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Figure 6. Correction for the continuum over-subtraction in the case of the 
cluster line emitter HI stack. 


HI gas than field line emitters (although not statistically significant) 
and about 12 times less than the cluster star-forming galaxies and 
AGNs (2.8(7 significance, where (7 is calculated as the errors on 
the passive galaxy mass and the cluster line emitter mass, added 
in quadrature). Note, however, that the velocity range used for the 
integration of the HI signal is larger for the passive cluster mem¬ 
bers than what was used for the line-emitters. The cluster emission¬ 
line systems have on average much lower masses compared to the 
passive galaxies. Therefore, the ratio of HI to stellar mass for the 
passive cluster galaxies is less than 0.008, a factor of~ 30 times 
lower than the cluster line emitters. Taking into account the errors 
on the gas fractions for the two populations summed in quadrature, 
the difference between the gas fraction in passive and cluster active 
galaxies is significant at the 4(7 level. 

3.6 Ha - HI correlation 

We compare the amount of cold gas and ionised content in each 
galaxy stack by investigating their Ha (Lh^) and HI (Lhi) lumi¬ 
nosities. 

The HI luminosity is calculated from the peak flux. The Ha 
luminosity is calculated from the Ha flux estimated from the NB 
observations (see § |2.2| ), after correcting for extinction by Galactic 
dust, as well as for 1 mag for intrinsic dust attenuation ( [Sobral 
|2012| ) within each galaxy. We also remove the contribution of the 
adjacent [Nil] line from the line flux (for details please see |Stroe| 
|et al.|2014^|2015| ). We average the corrected Ha luminosities for 
galaxies within each stack to obtain a mean value for the ionised 
gas content as function of galaxy type. 

The luminosities are calculated in the following way: 

^Hl,Ha = 47rJL^Hl,Ha (3) 

where F is the Ha total flux and HI peak flux, respectively and 
dh = 940 Mpc is the luminosity distance at the redshift of the clus¬ 
ter. 

The Ha can be converted into a SFR value using the rela¬ 
tionship from |Kennicutt| ( 1998| l, which we correct for the |Chabrier| 
( |2003| ) IMF, according to Salirn et ar] ( |2007| ): 

SFR ^ 4.4 X IQ-^^Lhc 
M0yr“l ergs“^ 

Note that not all HI stacked galaxies have an Ha flux measure¬ 
ment (see Table 1^ for numbers). This may be because the sources 
have too faint Ha line fluxes, below the limits of our NB Ha sur¬ 
vey. To test how the full HI sample differs from the HI sample with 
Ha measurements, we followed the procedure outlined in §3.2|3.4| 
and stacked only the HI sources with Ha measurements. We find 
that the peak HI fluxes and the average HI masses for subsamples 
with Ha counterparts matches their parent sample within the error. 


Table 5. Optical properties of each galaxy stack. The average Ha luminosi¬ 
ties for each stack are calculated from data from |Stroe et al.|t2015) . The 
SFRs are calculated using the (Kennicutt|1998) c onversion. Stellar masses 
are obtained using the method from |Sobral et al.|{2015) . 


Sample 

Number 

logio(THa) 
(erg s“i) 

SFR 

(M© yr-l) 

M* 

(10® Mq) 

Emission line, field 

22 

41.53 

1.49 

4.8±0.8 

Emission line, cluster 

45 

41.45 

1.23 

7.4±0.5 

Passive, cluster 

154 

<40.9 

<0.35 

25.6±0.5 


Therefore, the subsamples with Ha measurements are representa¬ 
tive of the parent sample. Given the more robust measurement of 
the average HI properties for the full HI samples (driven by the 
higher number statistics), in comparing the Ha and HI properties, 
we use average HI properties derived for the full samples. 

As Figure]^ shows, the Ha line emitters tend to be more lu¬ 
minous in HI. This is equivalent to galaxies which are more star¬ 
forming possessing larger reservoirs of atomic gas (Figurej^. 

3.7 Ha - radio correlation 

We extract sources from the VLA 1.5 GHz image using PyBDSM at 
the positions of the passive and Ha line emitter galaxies with spec¬ 
tra (for number of sources see Tablej^. The software detects single 
sources as islands and fits the flux distribution with Gaussians and 
also calculates the background noise levels using emission-free re¬ 
gions of the sky nearby each source. We assign a source a radio flux 
density by summing up the flux from all Gaussians belonging to its 
island. We cross-match radio sources with optical counterparts in 
our optical spectroscopic catalogue, using a maximum search ra¬ 
dius of 5 arcsec, to account for the positional accuracy of the op¬ 
tical and radio images as well as any extent the radio sources may 
have. 

In case a source is not detected in the radio map, we assign it 
an upper limit flux equivalent to 3(7 rms^ where the cjrms is calcu¬ 
lated from the noise level at the position of the source. Note that the 
FOV of the VLA image is large enough (FWHM of ~ 30 arcmin) 
that is covers all the optical source positions. 

We calculate observed 1.4 GHz measurements from the 1.5 
GHz values assuming a —0.7 radio spectral index value, and then 
convert the values to restframe 1.4 GHz measurements. 

A plot showing the relationship between the Ha luminosi¬ 
ties and 1.5 GHz luminosities (calculated using equation can 
be found in Figure The fluxes of the radio counterparts and 
their morphologies can be found in the Appendix in Table ] A l] The 
emission-line galaxies have, on average, 1.5 GHz luminosities 1—2 
orders of magnitude lower than the passive galaxies. Interestingly, 
even though the emission-line cluster members have similar Ha 
luminosities, and hence SFRs, to the field line emitters, their ra¬ 
dio BB detection rate is a factor of > 5 higher. However, for all 
emission-line sources with radio detections, the Ha luminosity cor¬ 
relates with the amount of radio emission. 


4 DISCUSSION 


Despite being extremely massive (M 200 ^ 2 x lO^^M©, 

Jee et al. 

2015| Dawson et al.||2015|) and hot (T = 6 — 12 keV, |Ogrean et 

al.|2013|), the ‘Sausage’ merging cluster hosts numerous 

massive. 


Ha emission-line galaxies displaying elevated levels of SF and 
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log(LHi) (erg s-^ Hz"^) 

Figure 7. Relationship between Ha luminosity and peak HI luminosity. 
As expected the cluster passive galaxies, selected to be non-line emitters 
contain less HI than emission-line galaxies. Line emitters are selected to be 
Ha luminous, indicating the presence of vigorous SF and/or AGN activity. 

AGN activity, outflows, high metallicities and low electron den¬ 
sities compared to galaxies in the field ( |Sobral et al.|2015[|Stroe et| 
|al.|20I5| ). In the present study, we find that the emission-line cluster 
galaxies have similar HI neutral gas fraction as the field galaxies. 
The data reveal linear correlations of the HI, Ha line emission and 
radio BB luminosity (Figures [^|^[^. By combining tracers of SF 
on different time scales, HI, Ha and radio BB data, we can under¬ 
stand the circumstances under which the elevated activity can be 
triggered and also the possible future evolution of the SF properties 
in the cluster galaxies. 


4.1 HI & Ha - tracing the gas that fuels future star 
formation episodes 

We make a clear detection of HI for the emission-line cluster galax¬ 
ies giving an average mass of (2.50 ± 0.62) x 10^ M0, while the 
average HI mass for the field galaxies is (1.86 ± 1.20) x 10^ M©. 
[Stroe et ^ ( |2015| ) and jSobral et al.j ( |2015| ) find that the stellar 
masses of Ha cluster galaxies are on average higher than their 
field counterparts. For the samples used in the HI stacks, the clus¬ 
ter galaxies are about 1.5 times more massive than the field line 
emitters (see Table [^. 

Note that the cluster and field line emitters are selected in the 
same way and that the spectroscopic samples are representative of 
their parent samples (see jSobral et al.|2015| ). Therefore the cluster 
and field emission-line galaxy samples within the ‘Sausage’ field 
are fully comparable. 

As mentioned in §[^ previous studies of the HI content in clus¬ 
ter galaxies find that star-forming galaxies become increasingly HI 
deficient towards cluster cores, when controlling for stellar mass or 
optical disk size (e.g. jSolanes et aL|2001[|Verheijen et aL|2007[|Lah| 
jet al.|2007] l. Contrary to previous work in the field, we find that our 
emission-line cluster galaxies are as gas rich as their field counter¬ 
parts, despite the two samples having similar Ha luminosities and 
hence similar SFR (see Table ^and F igu re]^. 

The studies of [Lah et ^ ( [2009| ) and jL^ et al.j ( [2007] ) indicate 
that a cluster at z ~ 0.37 and a blank field at z ~ 0.24 follow a sim¬ 
ilar relationship between SFR and HI mass to local field galaxies 
( [Doyle & Drinkwater|2006| ). Even though our galaxies show evi¬ 
dence for a correlation between the amount of Ha emission (or the 
SFR) and the HI mass, both the emission-line and passive galax¬ 


' ■ ' Passive, cluster 

Emission line, field 


Abell 370, z = 0.37, Lah et al. 2009 
Field at z = 0.24, Lah et al. 2007 


' ® ' Emission line, cluster Field, z ~ 0.0, Doyle & Drinkwater 2006 



Figure 8. Relationship between SFR and integrated HI mass. All SFRs are 
calculated with a Chabrier IMF. For the passive population, the SFR is an 
upper limit. Overplotted are the data points for cluster Abell 370 at z = 0.37 
jLah et al.|2009) and a blank field at z = 0.24 jLah et al.|2007t, together 
with the SFR-M/// relationship for the local Universe ( [Doyle & Drinkwate^ 
[2006) . Note that it is difficult to compare our dataset with other work be¬ 
cause of the different ways of measuring the HI mass and the different SF 
tracers used (Ha, [Oil] or IR). 


ies do not follow the [Doyle & Drinkwater[ ( |200^ relationship (see 
Figure [^. The HI masses of our sample are > 5(7 away from the 
masses predicted by the relationship at the same SFR. This could 
be entirely driven by the different SF tracers used in the different 
studies ( [Doyle & Drinkwat^ ([20061) use infrared data, [Lah et al^ 
( [2009) use [Oil] and [Lah et al.[ ( [2007] ) the Ha emission line) or the 
spatial or velocity range over which the HI signal was integrated 
over. 

Given the massive cluster galaxies may reside in massive dark 
matter haloes, they could have retained their cold HI gas more eas¬ 
ily during the cluster merger. Interestingly, while spiral galaxies 
in the Virgo cluster are highly HI deficient, they are not deficient 
in molecular gas { [Kenney & Young|1986[[ 19891 [Stark et al.|1986[ l. 
The authors attribute this to the preferential stripping of low-density 
gas located at the galaxy outskirts, therefore not affecting dense 
molecular gas located towards the galaxy centre. Therefore, in the 
case of the ‘Sausage’ line emitters, with little to no ram pressure 
stripping of neutral and molecular gas, the larger reservoirs could 
fuel increased SFR in the cluster galaxies. If the cluster galaxies 
maintain their current average level of SF (^ 1.23 Mq yr“^, see 
Table [^, and assuming 100% efficiency in converting cold gas into 
stars, the HI reservoir would be depleted in ~ 2.0 Gyr. If we as¬ 
sume a molecular gas content equal to the HI mass, it would take 
about ^4.0 Gyr to consume the gas. However, as shown in|Sobral| 
[eTaTl l IM?] ), the cluster galaxies also lose gas through outfiows. 
Assuming a maximal mass outflow rate similar to the SFR ( [Forster[ 
[Schreiber et al.[[200^ [Genzel et al.[[201l| , the HI gas will have 
been used up in about ~ 1.0 Gyr (~ 2.0 Gyr if we include molecu¬ 
lar gas). Assuming a more realistic outfiow rate of about 0.1 — 0.5 
SFR (as observed by Swinbank et al. submitted), the HI gas would 
be depleted in 1.35 — 1.85 Gyr, or 2.7 — 3.7 Gyr if molecular gas is 
considered. This is in line with calculations from [Stroe et al.[ { |2015[ ) 
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Table 6. As Table|^but only for sources with Ha measurements. Optical and broad-band radio properties of the sources in the line emitter and passive stacks. 
The average Ha luminosities for each stack are calculated from data from |Stroe et al.|{T015) . The SFRs are calculated using the jKennicuttH 1 99^ conversion, 
with a |Chabrie^j2003) IMF. 


Sample 

Number HI sources with 
Ha counterparts 

(erg s“i) 

SFR' 

Mq yr-‘ 

Number Ha 

sources 

Number Ha sources with 
radio counterparts 

logio(i'Ha)^ 

(ergs”*) 

Emission line, field 

20 (~91%) 

41.53 

1.49 

39 

3 ( ~ 8%) 

41.30 

Emission line, cluster 

29 (- 65%) 

41.45 

1.23 

54 

24 (~ 44%) 

41.20 

Passive, cluster 

79 (~51%) 

<40.9 

<0.35 

90 

11 (- 12%) 

40.15 


^ Average over the HI stacked sources with Ha counterparts. ^ Average over all spectroscopic sources with Ha measurements. 


^ SF 
AGN 

Emission line, cluster 


Emission line, field 
' Passive, cluster 


Field at z = 0.24, Lah et al. 2007 
Field at z ~ 0.0, Sullivan et al. 2001 



Figure 9. Relationship between Ha luminosity and restframe 1.4 GHz ra¬ 
dio luminosity. The vertical black dashed line marks the RMS noise at the 
centre of the radio image, while the gray line represents the detection limit 
towards the edges of the FOV. The horizontal, solid, black line indicates the 
detection limit of the Ha NB survey. Passive galaxies, undetected in Ha, 
are only plotted if they have a detection in the radio. Error bars are plotted 
but in most cases they are smaller than the symbol size. On average, the 
galaxies undergoing SF and optical AGN episodes have lower radio lumi¬ 
nosities, probably driven by SNRs, while the passive cluster galaxies host 
powerful jetted and tailed radio sources. It is noteworthy that, in the case 
of line emission sources with radio detections the Ha luminosity correlates 
well with the radio luminosity. We overplot for comparison data points for 
held galaxies at z = 0.24 from |Lah et al.|{2007) and the linear correlation 
between Ha and radio luminosity for local galaxies found by|Sullivan et al.| 
( | 20 ^ . 

where the molecular gas content was estimated using the total stel¬ 
lar mass, but atomic gas was not taken into account. 


4.2 Radio broad data - tracing the SN emission 

As shown by |Sullivan et al.| ( |2001| ), synchrotron radio emission in 
star-forming galaxies is generated in super-nova remnants (SNR). 
Given the time required for a 7 — 8 M© star to evolve to the red 
giant phase and undergo core-collapse, SNRs are good tracers of 
SF episodes happening ~ 100 Myr ago { |Condon|1992| l. 

In the case of Ha luminous cluster galaxies, undergoing 


strong SF and optical AGN activity, the Ha luminosity correlates 
well with the radio BB continuum (see Figure|^and Table |Al| in the 
Appendix). The values fall on the same correlation as the large sam¬ 
ple (^ 350) of typical z = 0.24, Ha luminous field galaxies studied 
by|Lah et al. (|2007| t and follow the tight relationship found by |Sul-| 
|livan et aLl ~i 2001| ) for local field galaxies. For both samples, the 
galaxies classified as purely star-forming or as hosting an optical 
AGNs follow the Ha - radio correlation. However, the AGN dom¬ 
inated cluster galaxies are expected to posses reasonable amounts 
of SF, as indicated by the spiral arm patters in many of their hosts 
(for images see |Sobral et aL|20T5] l. Therefore, the sample includes 
photoionised broad line and narrow line regions ( |Sobral etal.|2015| ) 
hosted by spiral galaxies, constituting examples of Seyfert galaxies. 
Even though some of the cluster galaxies are currently dominated 
by AGN, they have undergone significant SF activity in the past. 

Despite their similar average Ha luminosities, the fraction of 
Ha luminous cluster members with radio counterparts is a factor 
> 5 higher than their field counterparts, down to a similar sensitiv¬ 
ity limit (see Tabl ean d Figure]^. This is consistent with results 
from |Sobral et al.| ( [2015| ) where they found evidence for strong out- 
fiows, probably driven by SN, only in the cluster galaxies and not 
in their field counterparts. Therefore, cluster galaxies have been un¬ 
dergoing SF for at least 100 Myr, while field galaxies have been rel¬ 
atively inactive or less active in the past. Increased SF episodes tens 
to a few hundred Myr ago triggered by the cluster merger and its 
associated shock would also lead to higher SN rates. Field galaxies 
have not undergone any interaction with the shock front, hence in 
their case, there was no trigger for SN. 

Overall, line emission galaxies both inside and outside the 
‘Sausage’ cluster follow the local relationship between Ha and ra¬ 
dio emission. Field emission line galaxies from a larger sample at 
z ~ 0.24 (~ 150 field galaxies [Lah et aL|2007| l fall on the same rela¬ 
tion. This indicates that the relationship between Ha and BB radio 
emission does not evolve from z ~ 0.2 to the present and that it does 
not depend on environment. The stellar populations in all galaxies 
that have been undergoing SF for longer periods ('^ 100 Myr), irre¬ 
spective of redshift or environment, seem to evolve similarly from 
the massive, short-lived stars which are responsible for producing 
the bulk of Ha emission to the slightly less massive stars whose 
explosions dominate the SN population. These results indicate that 
the star formation history for the ‘Sausage’ cluster galaxies is rela¬ 
tively constant, without any strong recent (~ 10 Myr) busts of star 
formation or in the past 100 Myr. 

A little bit over 10 per cent of the passive galaxies have a ra¬ 
dio BB counterpart, a similar rate to field emission line galaxies, 
much 4 times lower than the cluster line emitters. The passive clus¬ 
ter members that have radio counterparts are giant ellipticals host¬ 
ing radio jets and tails as indicated by the radio morphologies and 
luminosities. By contrast to the emission-line galaxies, where the 
radio emission is most likely produced by SNR, the radio emission 
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in the elliptical galaxies traces shock-accelerated electrons in the 
jets and their back-flow (see also the Appendix). 


4.3 Relationship to cluster merger state and shocks 


As the radio BB and Ha data indicate, the ‘Sausage’ cluster galax¬ 
ies have been undergoing intense SF and AGN activity for at least 
100 Myr and this is likely to last a further ~ 1 Gyr, given the large 
reservoirs of neutral hydrogen. The SF can last for an another ~ 2 
Gyr if comparable amounts of molecular gas are present. How¬ 
ever, in the held galaxies around the cluster, we only And evidence 
of very recent SF episodes 10 Myr). Despite the comparable 
amounts of HI, the cluster galaxies therefore underwent a signifi¬ 
cant event triggering SF about ~ 100 Myr ago, evolving differently 
than their held counterparts. The most significant difference be¬ 
tween the cluster and field line emitters is the cluster merger and 
the passage of the merger-induced shock waves only affected the 
cluster members. 

The ‘Sausage’ cluster is a result of a massive merger about 0.5 
Gyr ago (e.g. |van Weeren et al.|201 l| l, which pr oduced shocks tray - 
elling through the ICM at about 2500 km s“^ ( Stroe et al.|2014c l. 
As the HI data indicate, the massive cluster members seem to have 
retained most of their neutral gas during the merger. Given their 
travelling speed, we expect the shocks to have traversed most clus¬ 
ter galaxies about 100 — 300 Myr ago. The SF time scale imposed 
by the radio and optical tracers therefore matches well with the 
cluster merger time line. Our results fully support the interpretation 
previously proposed by |Stroe et al.| ( |2015] ) and |Sobral et al.| ( [T015| ), 
where the cluster merger induced shocks trigger gas collapse as 
they traverse the gas rich cluster galaxies. This interpretation is also 
supported by simulations ( |Roediger et al.|2014] |. 


5 CONCLUSIONS 

We presented deep HI observations combined with Ha and broad 
band radio data to study the past, present and future SF activity in 
the ‘Sausage’ merging cluster. Our main results are: 

• The cluster Ha emission-line galaxies (star-forming and 
radio-quiet broad and narrow line AGN), selected down to the 
same SFR limit, have as much HI gas ((2.50 ±0.62) x lO^M©) as 
the held counterparts around the cluster ((1.84 ± 1.20) x lO^M©), 
when accounting for the different stellar masses of the two sam¬ 
ples. This indicates the massive cluster line emitters retained their 
gas during the cluster merger. 

• A stringent upper limit is placed on the average HI content of 
the passive galaxies in the ‘Sausage’ cluster: Mhi = (0.21 ±0.15) x 
IO^Mq. The ratio of HI to stellar mass for the passive galaxies is 
almost 40 times less than for cluster line emitters (significant at 4(7 
level). 

• If the present SF and outflow rate is maintained in the 
emission-line cluster galaxies, their HI reservoirs will be depleted 
in-0.75-1.0 Gyr. 

• A large fraction of the emission-line cluster galaxies have ra¬ 
dio BB detections, indicating the presence of SNR. These sources 
have been therefore undergoing vigorous SF for at least 100 Myr. 

• The relationship between Ha and radio continuum emission 
shows no evolution from z — 0.2 to the present and also no depen¬ 
dence on environment. 

Our HI observations represent an important milestone in the 
study of the ‘Sausage’ cluster SF history. The member galaxies are 


gas-rich (gas to stellar mass ratio of — 0.34) and thus capable of 
sustaining the increased SF and AGN activity measured in the clus¬ 
ter. 
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Table Al. Table with the radio fluxes, errors and morphologies of the VLA 
1.5 GHz counterparts to the spectroscopic sources. The same sources tabu¬ 
lated here are plotted in Figure]^ 


RA 

DEC 

Flux 

Error 

Morphology 

(deg) 

(deg) 

(mJy) 

(mJy) 



Passive, cluster 


340.6978 

53.0939 

15.034 

0.022 

NAT 

340.7687 

53.1234 

0.023 

0.006 

unresolved 

340.6048 

52.9719 

6.304 

0.033 

WAT 

340.7073 

53.0401 

0.034 

0.007 

unresolved 

340.8293 

53.1228 

12.420 

0.094 

WAT 

340.7072 

53.0081 

5.109 

0.125 

unresolved 

340.7051 

53.0920 

0.160 

0.013 

unresolved 

340.7191 

53.0806 

90.106 

0.090 

tailed 

340.7132 

53.0139 

8.084 

0.053 

WAT 

340.7385 

53.1344 

0.124 

0.008 

unresolved 

340.7720 

52.9552 

0.103 

0.010 

unresolved 


Emission line, held 


340.3917 

53.0163 

0.194 

0.014 

unresolved 

340.6917 

52.8561 

0.404 

0.025 

unresolved 

340.9183 

52.8727 

0.097 

0.013 

unresolved 


Emission line, cluster 


340.4575 

53.0954 

0.138 

0.011 

unresolved 

340.5930 

53.0083 

0.029 

0.006 

unresolved 

340.6151 

52.9683 

0.185 

0.011 

unresolved 

340.6518 

53.0904 

0.073 

0.007 

unresolved 

340.6598 

53.0138 

0.100 

0.010 

unresolved 

340.6642 

52.9674 

0.361 

0.012 

unresolved 

340.6711 

52.9746 

0.817 

0.026 

disk 

340.7136 

52.9061 

0.133 

0.015 

disk 

340.7225 

52.9517 

0.080 

0.009 

unresolved 

340.7369 

52.9439 

0.439 

0.015 

disk 

340.7401 

53.0874 

0.336 

0.022 

disk 

340.7475 

53.0889 

0.179 

0.017 

disk 

340.7502 

53.0439 

0.139 

0.010 

unresolved 

340.7796 

52.9257 

0.122 

0.013 

disk 

340.7870 

53.0902 

1.207 

0.018 

unresolved 

340.7941 

53.0716 

0.135 

0.011 

disk 

340.7985 

53.0764 

0.055 

0.008 

unresolved 

340.8037 

53.0029 

0.285 

0.016 

disk 

340.9204 

53.0778 

0.344 

0.016 

disk 


NAT: narrow angle tailed galaxy; WAT: wide angle tailed galaxy; 
Disk: galaxies with a disk or spiral-like morphology. 


APPENDIX A: RADIO FLUXES AND MORPHOLOGIES 

We tabulate the flux values of the 1.5 GHz VLA counterparts to the 
spectroscopic sources in Table |Al| We also describe the morphol¬ 
ogy of the radio sources. Passive galaxies are hosts to jetted radio 
AGN, mainly pushed in wide or narrow angle tail morphologies by 
the interaction with the ICM. Emission line galaxies have mostly 
disk or spiral-like morphologies, indicating the radio emission is 
coming star formation. 








